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The main strategies to reduce lattice artifacts for spin models, gauge fields, free fermions and QCD are discussed. 
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1. INTRODUCTION 

The idea to improve lattice actions is 
nearly as old as the use of lattices in regularizing 
field theories j^]. Nevertheless, until recently nu- 
merical simulations were performed mostly with 
the simplest (standard) form of lattice actions. 
This attitude seems to have been altered — at 
this conference there were probably more contri- 
butions using some kind of improved action than 
the standard ones. The reason is, obviously, that 
the lattice artifacts resulting from the finite lat- 
tice spacing a cannot be easily reduced further 
simply by brute force, i.e. increasing the lattice 
size. As explained in Ref. Q the computing cost 
for full QCD grows at least as I would like 

to add to this that in some cases the situation is 
much worse, even for pure gauge theory. 

• Thermodynamics. 

For measuring the basic thermodynamic 
quantities, like the free energy, a conserva- 
tive estimate gives 



computing cost 
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(This is because the free energy density, 
/ = 0(T'*), is hidden behind the uhra- 
violate fluctuations of the action density, 
(S) = 0(l/a^), as will be discussed later.) 

• Topological susceptibility. 

Using the standard form of the action (with 
the geometrical definition of the topological 
charge) the presence of dislocations — con- 
figurations with topological charge and too 
small value of the action - spoil the result 
even in the continuum limit a ^ 0! 



In this talk I will concentrate mostly on theo- 
retical ideas and developments, since the numer- 
ous contributions concerning simulations with im- 
proved actions are covered in several plenary talks 
at this conference. 

2. IMPROVING THE ACTIONS 

There are two main ways of improving the lat- 
tice actions, based on: 

• Wilson's renormalization group (RG), 

• Symanzik's approach, inspired by perturba- 
tion theory (PT). 

Wilson's approach |^ is based on the fact that a 
critical system (with a correlation length ^ = oo) 
remains critical after mapping the lattice to a 
coarser one by a block transformation. The lat- 
tice action itself changes, in general, under the 
renormalization group transformation (RGT) and 
moves (except for pathological transformations as 
decimation) towards a fixed point (FP) on the 
critical surface. The trajectory leaving the crit- 
ical surface along the relevant (unstable) direc- 
tion(s) defines the renormalized trajectory. The 
system specified by an action on the renormalized 
trajectory has exactly the same partition function 
as the original system defined on an 'infinitely 
fine' lattice, and describes the same long-range 
physics with no lattice artifacts. The RGT can 
be (to a large extent) arbitrary, leading to differ- 
ent FP's and renormalized trajectories. There is, 
however, an additional useful option — one can 
vary (and optimize) the RGT during the succes- 
sive RG steps on the way to smaller and smaller 
correlation length. One can call this a 'running' 
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RGT. The actions in this large set have the prop- 
erty that they are perfect - i.e. produce no cut-ofF 
effects (for physical quantities like masses, scat- 
tering amplitudes). 

The disadvantage of this beautiful construc- 
tion is that the actions involved have complicated 
structure, and it is difficult to determine and use 
them in numerical simulations. Perhaps to avoid 
these difhculties, Symanzik |^ suggested an alter- 
native approach. The procedure was to add some 
extra terms to the naive lattice action and choose 
their coefficients appropriately to cancel the lead- 
ing O(a^) lattice artifacts in the n-point Green- 
functions up to a given order in PT. (O(a^) refers 
to the bosonic case; for Wilson fermions the lead- 
ing artifacts are more severe, 0(a)). For gauge 
theories this requirement had to be reduced to 
'on-shell improvement', i.e. to the perturbative 
improvement of physical quantities The 
number of extra terms to be added is determined 
by the number of possible dim = (5), 6 operators. 
By allowing redefinitions of fields, however, this 
number can be reduced. The number of extra 
operators to be added is following: 



op s 



order Ref. 



gauge fields: 
fermions 
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(Some of the 10 fermionic operators contributing 
to O(a^) are quartic in the fermion fields.) The 
statement is that to arbitrary order in only 
the coefficients of these operators have to be de- 
termined to cancel the corresponding cut-off de- 
pendence of 0((7"a) or 0(5"a2). Although the 
Symanzik approach provides a relatively simple 
modification of the action to cancel the leading ar- 
tifacts, the perturbatively determined coefficients 
are usually not applicable in the case where most 
of the numerical simulations are done. The other 
point, to which I shall return, is that there is a 
large freedom of choosing the given number of 
extra operators, and some additional criteria are 
needed to find the choice best suited for numerical 
simulations. 



3. RECENT DEVELOPMENTS 

I will consider three main recent developments 
in the improvement program: 

• non-perturbative numerical determination 
of 0(a) improvement coefficients in QCD 

i 

• tadpole improvement (TI) |9|Jl0t] 

• FP actions (or classically perfect actions) 



The first one is a direct application of 
Symanzik's improvement program for the non- 
perturbative situation in QCD: instead of cal- 
culating the 0(a) cut-off effects for some quan- 
tities in PT, the authors measure these effects, 
and by tuning the coefficient csw idn) of the 
Sheikholeslami-Wohlert (clover) term [Q they get 
rid of the leading 0(a) artifacts. Since it is done 
in the framework of the full QCD, I will discuss 
it at the end, and start with the simpler case of 
bosonic theories. 

4. TADPOLE IMPROVEMENT 

The motivation behind this approach is the ob- 
servation that large coefficients appearing in PT 
come dominantly from tadpole diagrams. Lepage 
and Mackenzie suggested ||] that these contribu- 
tions should be removed by a mean-field redefini- 
tion of the original fields. The corresponding pre- 
scription for the SU(3) gauge fields is to replace 
the fields in the original lattice action according 
to the rule 

U^{x) ^ -U,,{x) , (1) 



Uq 



where 



(2) 



In addition, to get a better perturbative expan- 
sion, it is suggested to use 

= -3:0^39 ^°sQ(C/pla,)) (3) 

as an expansion parameter instead of aiatt = 
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For tadpole improving (TI) a one-loop 
Symanzik improved (SI) action (where the coef- 
ficients are g dependent) the prescription is to 
recalculate the improvement coefficients for the 
modified action obtained by (||) using ordinary 
PT, substitute aiatt ois, and use this action in 
the numerical situation. The procedure involves a 
self-consistent determination of the coefficients in 
the action: the coefficients to be used in the MC 
simulation have to be determined from MC simu- 
lations. Technically, this is not a serious problem, 
since it is very easy to measure (?7piaq). (Princi- 
pally, however, the mean field definition can cause 
some problems. For example, (C/piaq) depends 
not only on (3 but also on the size of the system. 
Should one change mo with the temperature at a 
given (3 according to eq. (||) or should it be deter- 
mined from the zero temperature, infinite volume 
case? One can argue both ways.) 

I cite here two examples of tadpole improve- 
ment. For the 0(a) one-loop Symanzik-improved 
Wilson fermions the coefhcient IIISNitI 



csiv = 1 + 0.2659g2 



(4) 



should be compared with the TI tree-level 
Symanzik-improved action which gives 



= ( 3RcTr(C/piaq) 



1 



-3/4 



= 1 



In other words, 94% of the original one-loop 
correction is removed by making the TI in the 
tree-level formula. Accordingly, tadpole improv- 
ing the 1-loop formula amounts to using Cgy/ = 
1 + Att ■ O.OlSOa^. For the pure gauge field TI 
for the 1-loop SI action of Liischer and Weisz 
(who have used the plaquette (pi) and the rect- 
angle (rt) ) means 



/3pi 



20 



(l-h2.0146a^) 



1 



(l+0.4805a. 



i.e. 76% of the coefficient has been removed by 
applying (|^) to the tree-level result. 

The TI action for gauge fields has been demon- 
strated to work well even for lattices as coarse as 
a ^ 0.4fm [Q. This success and the simplicity 
of the prescription resulted in a large number of 



contributions to this conference which apply TI 
actions in numerical simulations. However, as we 
shall discuss later, TI does not work everywhere, 
e.g. for full QCD it does not give properly the 
value of csw at the /? values of interest |^ . 

After discussing the FP actions, I will try to 
suggest how to understand the tadpole improve- 
ment from the corresponding point of view. 

5. FIXED POINT ACTIONS 

To be specific, I will consider here the pure 
gauge theory. A RGT is defined by specifying 
the kernel of the block transformation, T{V^U). 
Here V is a gauge configuration on the coarse lat- 
tice while U denotes a configuration on the fine 
lattice. Given the action S{U) one defines S'{V) 

by 



P'S'(V) ^ / ^^^-t3[S{U)+nVM)] 



(5) 



{S'{V), as well as S{U), is normalized to the 
continuum expression for very smooth fields.) 
T{V, U) should satisfy a normalization condition 
to leave the partition function unchanged. Since 
the theory is asymptotically free, for /3 ^ oo one 
has /3' — (3+0{l), and in this limit eq. (||) reduces 
to a saddle point equation: 



S'{V) = uAn{S{U)+T{y,U)} 



(6) 



The fixed point of this transformation is given by 



S 



FP 



[V) = mm {S'PP([/) -I- T{V, U)] 



(7) 



Note, that although the (3 ^ oo limit has been 
taken, the fields V and U are not assumed to 
be necessarily smooth — they could be arbitrar- 
ily rough as well in this definition. Eq. (0) can 
be solved analytically for smooth configurations 
and numerically for arbitrary V by iteration. The 
iterative procedure converges very rapidly since 
the typical action density on the fine lattice is 
^ 30 times smaller than that on the coarse lat- 
tice (a factor 16 coming from the volume, and 
typically half of the value of the r.h.s. is due to 
T{V,U).) Consequently, in most cases one can 
use on the r.h.s. just the quadratic approxima- 
tion to S^^ {U) or even the standard action. 
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The FP action has remarkable properties — it 
defines a classically perfect lattice action in 
the following sense. The solutions to the lattice 
equations of motion are related to their contin- 
uum counterparts, with the same value of the ac- 
tion: 



'S'cont(C4ont) — 'S'latt(^latt) • 



(8) 



In particular, for continuum theories possessing 
instantons, the FP action has scale invariant in- 
stanton solutions as well. One can say that the 
FP action is 'tree-level Symanzik-improved to 
all orders in a'. By proper choice of T{V,U) 
(which is achieved by tuning some free parame- 
ters in it) S^^ can be made quite short ranged — 
the couplings outside the unit hypercube can be 
suppressed (relative to the plaquette) by a factor 
^ 100 The FP action is not expected to be a 
perfect action — using it in numerical simulations 
with finite f3 (at finite correlation length) some 
cut-off dependence remains. Numerical simula- 
tions for 2d spin models and Ad SU(3) 
gauge theory show, however, that the lattice 
artifacts using FP actions are negligible for these 
theories even on extremely coarse lattices. 

Actions which are (quantum) perfect are ob- 
tained by following the renormalized trajectory 
(also allowing to vary T(V,U) with (3). Strong 
arguments suggest that the FP action is also 1- 
loop perfect, i.e. 



/^S-P^f ([/; (3) = pS^^iU) [1 + (/3-2)] 



(9) 



This conjecture is supported by formal RG argu- 
ments [p]Jll|] and by numerical evidence for the 
0(3) non-linear sigma-model |2^]. In |^l| it was 
found that the mass gap in a small box calcu- 
lated to 1-loop in PT with the FP action has no 
0((7^a") error. A coefficient R appearing in the 
1-loop expansion of the mass gap is plotted vs. 
(a/L)2 inFig.|l|. As seen, the correction is vanish- 
ing as exp(— 7L/a). This reflects the requirement 
that the size of the box should be larger than the 
interaction range of the FP action (otherwise the 
infinite-volume action should be modified). The 
tiny power-like artifact seen in the fit is caused by 
restricting the interaction distances considered in 
the perturbative calculation. 
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Figure 1. The value of R vs. (a/L)^ for the FP 
action for L/a = 2,3, ...,10. The fit shown is 
-0.1422022 - 0.0000176 (a/L)2. Note that the 
exact limiting value is —0.1422024. 



Although the FP action is classically and pre- 
sumably 1-loop perfect, the approximation actu- 
ally used in MC simulations contains some imper- 
fect steps: 

• truncation in the interaction distance, 

• parametrization (choice of the set of opera- 
tors to fit the numerically known values of 
the FP action) 

• the FP action is not on the renormalized 
trajectory (it is not 'quantum perfect'). 

One has to control these points to get a reliable 
FP action. It produces a misleading result, for ex- 
ample, to restrict the RGT to a small subspace of 
operators beforehand. One must show first that 
the interaction range is not larger than the trun- 
cation radius for the interaction terms and trun- 
cate the set of operators only afterwards. 

5.1. Natural choice of improvement terms 

In Symanzik's approach one has a large free- 
dom to choose the improvement terms. For the 
tree-level improvement to Oio?) one can choose 
the following loops, for example, 

A.) the plaquette and the rectangle (2x1 loop). 
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B.) the plaquette and the 2x2 loop. 

Which one should be preferred? Let me illustrate 
the point on the 2d tr-model. Denote by pxy the 
coupling between two spins at the relative dis- 
tance {x,y). (The standard action is given by 
the nearest-neighbour term pig = 1.) Take two 
choices of operators in the Symanzik approach: 

a. ) pio = 4/3, P20 = -1/12, 

b. ) pio = 2/3, pii = 1/6, 

The FP values for these coefficients (for an 'opti- 
mal' RGT 

pio = 0.618, pii = 0.190, P20 = 0.002 

The choice a.) has an antiferromagnetic coupling 
and a complex dispersion relation in the quadratic 
approximation. Clearly, our RGT prefers the 
choice b.) in the Symanzik program, which is 
free from these diseases. 

Similarly, for gauge theory the FP action in 
quadratic approximation gives a much smaller co- 
efficient for the rectangle (rt) than for the 'bent 
rectangle' (brt) or the parallelogram (pg) JTH . 
(Note that the brt and pg are equivalent in the 
quadratic approximation.) In the standard tree- 
level SI form, however, the plaquette (pi) and the 
rt are the two operators chosen to have non-zero 
coefficients. This choice does not correspond to 
the FP action. In quadratic approximation it pro- 
duces again a complex spectrum above some mo- 
menta. In fact, Liischer and Weisz Q proposed 
a one-parameter family of SI actions with a free 
coefficient x for the brt, suggesting that x ^ 
could be useful in MC simulations. However, as 
far as I know, only the x = choice has been 
exploited. 

I think, for the Symanzik program a natural 
choice of the improvement operators could be 
the one inspired by the KG approach. One can 
also turn the argument around. As mentioned 
above, the FP action is automatically Symanzik- 
improved, to all orders in a. A necessary trun- 
cation to a reasonable number of operators does 
not keep, of course, this property. Nevertheless, 
for a FP action with a short interaction range 
the violations of the lowest order (in a) Symanzik 



conditions will be small. Since the truncation 
is to some extent an arbitrary procedure, one 
might also keep these conditions satisfied. I note 
here, that we were interested mainly in a good 
parametrization of the FP action at relatively 
rough configurations (typical at short correlation 
length) hence we did not always impose the above 
conditions. 

5.2. Parametrization of the FP action 

In jl^] the FP action has been parametrized as 

^^^(c/) = EE^"(^)^", (10) 

C m 

with 

Xc =ReTr(l-C/c) , (11) 

and C denotes a closed loop. A few simple loops 
were taken with their powers m — 1..4. The 
higher powers are essential in this parametriza- 
tion, and they do not cause much overhead in 
numerical simulations. With a Swendsen-type 
blocking in ref. ||l^ a good scaling was observed 
for a < 0.4fm. In addition, a good asymptotic 
scaling for Tc/A has been found, which could in- 
dicate that the bare PT for the FP action works 
better. 

On the renormalized trajectory the coefficients 
in eq. ( |lO| ) will be g dependent. The 1-loop per- 
fectness, however, predicts that this dependence 
is given by 

Cm{C,g)^CmiC) + 0{g^), (12) 

i.e. no 0{g^) corrections are expected. One con- 
cludes that no large tadpole contributions could 
be present for the FP action, and it would be 
incorrect to tadpole improve it by 

^FP(f/)^^FP(C//yo). (13) 

On the other hand, a mean field approximation 
to the FP action ( |lO| ) yields an action of the TI 
type: 

5^P«E^<=ff(C'(^c))^c, (14) 
c 

where c,s{C, (Xc)) ^ ci{C) + 2c2{C){Xc) + ... In- 
deed, the relative weight of the parallelogram (pg) 
to the plaquette, Ccff (pg)/coft(pl) grows with g^ as 
suggested by the tadpole improvement (where it 
is oc 1/uq). 
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5.3. Improving correlators of local opera- 
tors; the qq potential 

In the Symanzik approach one improves the 
action to cancel cut-off effects in spectral quan- 
tities. To achieve this for correlation functions 
one also has to improve the operators. The sit- 
uation is the same in the RG approach. For a 
square-shaped blocking, although the spectrum 
is exact, the potential shows violations of rota- 
tional symmetry. (A close analogy: in the con- 
tinuum electrodynamics the interaction between 
square-shaped charges.) The perturbative qq lat- 
tice potential has, in general, terms of the follow- 
ing type: 



V{r) 



-P4(cOS( 



The second term (related to the distorted gluon 
spectrum) is removed by improving the action, 
the third one (an octupole term) by improving 
the operator. In the RG approach, the last term 
comes from the interaction of the quantum fluc- 
tuations on the fine lattice for a given coarse con- 
figuration (and is similar in nature to the expo- 
nentially vanishing term in Fig. In Ref. ||Tl[] 
using a Swendsen-type RGT (our type-1 RGT) 
we have got a qq potential which gave the proper 
string tension at a = 0.4fm but had consider- 
able violation of rotational symmetry at small dis- 
tances. This has been associated in ref. with 
the strongly square-shaped averaging of the RGT 
applied. To illustrate this point, and also to ob- 
tain an even better FP action, in Ref. |Q a more 
general RGT (type-3) has been considered, which 
included generalized staples along the 2d and 3c? 
diagonals. As expected, the rotational invariance 
at small distances has been improved both for the 
perturbative and the non-perturbative potential. 
In Fig. H the deviation of the perturbative po- 
tential from the continuum result is shown. In 
Fig. ^ the qq potential measured at finite tem- 
perature T = 2/3 • Tc and a « 0.4fm is com- 
pared for the Wilson action and type-3 FP action 
||2^ . For comparison, I measured also the 1-loop 
tadpole-improved potential for this temperature, 
also shown in Fig. ^. Surprizingly, the tadpole 
and the type-3 FP data almost coincide. 

Although it is not related to the FP actions, I 
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Figure 2. The difference V{r) — ycont('') vs. r for 
the Wilson action, type-1 and type-3 FP actions. 



end this section with a reference to where it 
has been suggested to test improvement schemes 
by using analytic and numerical methods in a 
small physical volume. The authors proposed a 
combination of the plaquette, the 2x1 and the 
2x2 loops, with coefficients 16/9, -1/9 and 1/144 
which also satisfy the tree-level Symanzik con- 
ditions. The improvement for the small volume 
quantities is about the same as for the Liischer- 
Weisz choice (5/3, —1/12 and 0). The choice of 
p3|] is motivated by the observation that their 
analytic calculations could be extended to inves- 
tigate the effects of the TI in this case. 

Let me also mention a contribution ^J] where 
the authors used the 1-loop SI TI gauge action 
on anisotropic lattices to calculate the glueball 
mass. The use of an anisotropic grid which is 
finer in the time direction (at < a^) is a useful 
tool to study heavy objects. The data for the 0++ 
and 2++ glueballs scale well up to Og — 0.4fm, 
while the standard action shows large artifacts 
on coarse lattices. 



7 



4.0 



2.0 



0.0 - 



-2.0 



-4.0 



T=2/3T, 


X 


xx^l* - 

X ± 




x> 

X 








X Wilson 
• type III 
O tadpole 



0.0 



0.5 



1.0 



1.5 



2.0 



Figure 3. V{r;T)/Tc vs. rT^ for the Wilson, 
type-3 FP and 1-loop TI actions, at finite tem- 
perature, T = 2/3 ■ Tc- The dashed line is the 
universal curve. 



6. THERMODYNAMICS 

Thermodynamic quantities measured at rela- 
tively large temperatures are sensitive not only to 
the low energy part of the spectrum hence the lat- 
tice artifacts in this case are usually much larger. 
For example, in PT for T ^ oo at temporal lat- 
tice size iVt = 4 the cut-off effect for p/T* using 
Wilson's gauge action is « 50% p5| . 

The pressure is given by 



rp4 



Nf / dp'{{S)^-{S)T) 



(15) 



l3o 



Here the action densities on the r.h.s. are mea- 
sured at a given /?' for two different temporal 
geometries, at A^t = oo (or rather for a 4-cube 
Nt = Ns) and Nt corresponding to the given 
temperature T. (Actually, this is the formula 
which shows that the corresponding computer 
cost grows like N}^ , as mentioned in the intro- 
duction.) In Ref. [g5| the authors measured p/T^ 
with the standard action for At — 4, 6, 8 and a 
few improved actions, including the 2x1 1-loop 
TI action for Nt = 4. (For the results see the 
figure in Ref. ||2^].) Papa has measured the 
same quantity for the FP actions with RGT of 



type-1 and type-3, for A* = 2,3 and T = 2Tc. 
The relative deviations from the continuum result 
(obtained by extrapolating the standard action 
results) for T = 2Tc are given below. 

At Wilson type-1 type-3 tadpole 



100% 



20% 



25% 
10% 



10% 

"0%" 



'0%" 



"0%" means that no deviation is seen within the 
error. Note that the At = 2 case is special be- 
cause of possibly large exp(— 7At) type correc- 
tions. 

One can conclude that the cut-off effects are 
drastically reduced using these improved actions. 

7. TOPOLOGICAL SUSCEPTIBILITY 

The FP action is especially suited for measur- 
ing topological effects: 

• it has scale invariant instanton solutions 
(for sizes p > 0.7a) 

• a good definition of topological charge is 
possible (classically perfect Q) 

Definition: Measure the topological charge on 
the minimized fine configuration U = U (V) of 
eq. (0), 



Q{V)=Q{U{V)) 



(16) 



(Here I mean a minimization iterated a few times 
if necessary, on even finer grids. Since on the 
fine grids the fields are rather smooth, after the 
second iteration the value of Q usually does not 
change.) The charge Q{V) avoids the problem of 
dislocations: any configuration V with Q{V) — 1 
has a value of the FP action not smaller than the 
continuum result, S^^ (V) > Sj^^V 

By measuring the topological charge this way, 
one also avoids problems associated with the cool- 
ing procedure (e.g. the disappearance of the in- 
stantons) — the minimization procedure provides 
a prescription to separate the fluctuations from a 
smooth background solution. The configuration 
U{V) is smoother than V, it is closer to a solu- 
tion — part of the energy of the fluctuations in 
V is taken over by T{V, U) in eq. (0). If F is a 
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solution then T{V, U{V)) — 0. Consequently, the 
minimization V —>■ U(V) provides an ideal cool- 
ing procedure. A serious drawback is that the 
minimization procedure is very slow, it is hard to 
use it in a MC measurement. A solution is to 
find an approximate analytical form for the func- 
tion U{V) as has been done for the spin models 

El®- 

This strategy has been applied to the 2d non- 



linear sigma- model yj, the 2d CP-" model |20| 



and to the 4d SU(2) gauge model |29[. For the 
0(3) model Xt^"^ does not scale (not unexpect- 
edly) while the CP'^ model shows a beautiful scal- 
ing, for ^ > 7. In addition, asymptotic scaling is 
found in m/A2_ioop for C ^ 4, suggesting again 
that the bare PT should work better for the FP 
action. Of course, in the 4d SU(2) gauge theory 
the technical difficulties are more serious, and the 
investigations are not completed yet. Neverthe- 
less, the results for the scaling of the topological 
susceptibility ||2^ are encouraging. For the corre- 
sponding plots see Ref.'s [ po|j2S|] . 

Let me discuss a few points in some detail. 
Blocking down a smooth instanton solution to a 
coarser lattice remains an instanton solution un- 
til p ^ 0.7a where p is the size of the instan- 
ton. Instantons of smaller size 'fall through the 
lattice' — they cannot be distinguished from lo- 
cal excitations. This is the consequence of the 
classical (saddle point) approximation. By fixed 
coarse configuration V at finite f3 in eq. (^) the 
typical fine configuration U fluctuate around the 
minimizing configuration U{V). These can have 
different charges, especially if the action density 
of V is localized to a small region. In contrast, 
U{V) assigns a unique fine configuration to V, 
and for instantons which 'fell through the lattice' 
this gives zero topological charge. Hence the con- 
tribution of these small instantons is lost by using 
the FP action and FP topological charge on the 
coarse lattice. Nevertheless, the definition (|l^) of 
the FP-charge associates Q = for these small 
size excitations so they don't spoil the topological 
susceptibility. One expects to observe scaling of 
Xt£,'^ once the role of instantons of size p < 0.7a 
becomes negligible. By using the standard ac- 
tion and even the (tree-level or 1-loop) SI action, 
there is a wide range of instanton sizes where the 



instanton action is smaller than the continuum 
value while the charge associated is 1 in the ge- 
ometrical definition — see the corresponding fig- 
ure in Ref.|Q. As a consequence, one can grossly 
overestimate xt and observe scaling only at much 
finer lattices. I also have the feeling that using 
TI will not help much in this case — the mean 
field approximation represents the action at some 
average 'roughness' of the fields, while instantons 
of different sizes are perhaps important at a given 
value of (3. In Ref. the authors discuss the no- 
tion of the 'quantum perfect topological charge' 
on the example of the exactly soluble Id XY- 
model. 

Note that the use a parametrized form for the 
minimized fine field U{V) technically might be 
similar to applying appropriately smeared fields 
to measure the topological charge pO[ |. 

8. NGN PERTURBATIVE DETERMI- 
NATION OF IMPROVEMENT COEF- 
FICIENTS 

In QCD with Wilson fermions the Wilson term 
(which has been introduced to avoid the dou- 
bling problem) explicitly breaks chiral invariance 
and induces 0(a) cut-off effects in chiral rela- 
tions. To remove these artifacts it has been pro- 
posed to modify the fermion action by a Pauli 
term acswi"^piu'4^Fiiu- The lattice representation 
of this operator is chosen to be the 'clover' term 
0. The coefficient csw(5o) has been determined 
in the SI approach to 1-loop order and given 
by eq. (|) (l|J|. In Ref. § the dependence 
of this coefficient on the coupling has been 
determined by non-perturbative methods. The 
framework used was the Schrodinger functional 
method, where instead of periodic boundary con- 
ditions in the time direction Dirichlet b.c.'s are 
used for the gauge fields and fermions. This has 
the advantage that a background field can be in- 
troduced in the system and the response to this 
field can be investigated. By using a small phys- 
ical box and Dirichlet b.c. the zero modes in 
the fermion matrix are avoided, hence this frame- 
work is especially suited to study the chiral limit 
m 0. 

From the PCAC relation d^A^^{x) = 2mP''{x) 
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for the axial current and pseudoscalar density, 
valid in the continuum limit, the authors define 



(17) 



where is a fermionic operator on the bound- 
ary. For Wilson fcrmions this quantity depends 
very strongly on the position x and the back- 
ground field, even at a = 0.05fm! Introducing 
the SW term into the action, and an extra term 
in the current via 



Al{x)^Al{x)+CAad^P''{x)., 



(18) 



by properly tuning the coefficients csw and ca 
one can make the quantity ( p7| ) independent of x 
and the background field. Fig. | shows the de- 
pendence of Csw on the coupling 



result a,t P — 6 fails by about 15% and the dis- 
crepancy seems to grow rapidly with go. 

In a work presented at this conference |3l| 
Gockeler et al. used the value csw predicted by 
the ALPHA collaboration Their results at 
/3 = 6 shown on an APE plot seem to extrap- 
olate well to the physical N/p ratio. R. Gupta 
noted that minimizing the source-dependence in 
his measurements prefers a value of csw(5o) con- 
sistent with the ALPHA result 

9. FERMIONS 

For fermions, compared to pure gauge fields, 
the cut-off effects are larger and the computa- 
tional cost grows faster with the lattice size, hence 
it is especially important to optimize the choice 
of the lattice action used. 
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Figure 4. Dependence of cgw on the coupling, 



We can conclude that the SW choice for the 
improvement term works well, but the coeffi- 
cient is strongly non-perturbative in the inter- 
esting region of couplings. Perhaps one can try 
other schemes where the higher order (or non- 
perturbative) corrections are smaller. Note also 
that the tadpole improvement of the 1-loop SI 



9.1. D234(|) action 

In their contribution to this conference, Alford, 
Klassen and Lepage suggested a modification 
of the 'D234' fermionic action which they 
call 'D234(|)'. It is constructed using higher or- 
der lattice derivatives along the axis. In addition 
the authors also introduce anisotropic lattices, by 
taking different lattice spacings along the spatial 
and time directions, as /at ~ 2 — 3. With this 
additional trick the dispersion relation for free 
fermions becomes real and reasonably good. (For 
the isotropic lattice the dispersion relation is com- 
plex above ap k, 1.) Note that the anisotropic 
choice, besides curing the spectrum has the ad- 
vantage that one can measure time correlations 
on a finer scale which can be important for heavy 
objects. The authors have plotted the 'velocity of 
light' c(p), defined via E{pf = Ei^d) + c(p)2p2, 
measured for different actions in the quenched ap- 
proximation for TT and p. While the TI of 1-loop 
SI gluon+fcrmion action gives a large deviation 
from 1, the 'D234(|)' action on an anisotropic 
lattice produces good agreement with c(p) = 1. 

I would like to point out that choosing an 
anisotropic lattice can be a useful tool in itself for 
dealing with heavy objects, but it is not needed to 
improve the spectrum of free fermions. As inves- 
tigations of the FP fermionic actions show this 
goal can be achieved by using appropriate cou- 
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plings on the unit hypercube. 

9.2. FP actions for fermions 

The RG approach to the fermion problem is 
presented at this conference by two groups, from 
MIT and from Bern/Boulder The ap- 

proaches of the two groups are very similar, the 
main difference is in the type of the blocking 
transformation used. This difference is of tech- 
nical nature, the results obtained and the expec- 
tations are very similar. 

The FP action for the massless free fermions 
has the form 

= E ^('') + Kr)] ^{n + r) . (19) 

The features of this FP actions are: 

• The spectrum is exact, there are no 0(a") 
cut-off effects ('classically perfect') 

• There are no doublers, but chiral symmetry 
is violated. This violation is, however, very 
special: it comes entirely from the blocking 
transformation which is not affecting the 
physics. 

• By optimizing the corresponding RGT the 
interaction range could be made quite small 
— couplings decrease exponentially fast and 
outside the unit hypercube could perhaps 
be neglected. 

The TO 7^ case is very similar, only in this 
case one does not have a FP, of coarse, since the 
mass is running. It is possible to start with an 
arbitrarily small mass and make an appropriate 
number of blocking transformations, doubling the 
mass in each step, until the given value of ma is 
reached. An important point here is the follow- 
ing: the RGT which was optimal for the to = 
case will not be optimal for the massive case — 
by applying a fixed RGT the interaction range 
would grow. To stay with a short ranged action 
one has to modify the parameters of the RGT 
with the growing mass. This can be done for all 
RGT's considered, and one can keep the action 
short ranged — at least up to ma « 2 but per- 
haps even higher. 



r 


A(r) 




0000 


2.312461 




1000 


-0.112098 


0.140881 


1100 


-0.034557 


0.034424 


1110 


-0.014912 


0.010742 


1111 


-0.007373 


0.003714 


2000 


0.001493 


0.001380 


3000 


-0.000088 


0.000243 


4000 


0.000000 


0.000010 



Table 1 

Some coefficients of a fermionic FP action. 



In numerical simulation one needs to trun- 
cate the interaction distance r in eq. (|9|). This 
introduces some imperfectness: the spectrum 
is distorted, 'ghosts' (higher lying and some- 
times complex branches in the spectrum) appear, 
Symanzik's improvement conditions for the spec- 
trum will be (slightly) violated. Nevertheless, 
since by the truncation one neglects only small 
coefficients, these effects are not expected to be 
large. Table |l] lists some of the coefficients of a FP 
action (corresponding to one of the Bern/Boulder 
RGT's) for free massless fermions. The coeffi- 
cients A(r) and pi{r) decrease rapidly with the 
distance r hence truncation to the unit hypercube 
seems to be justified. 

Fig. H shows the spectrum E{q, 0, 0) (along the 
axis) at ma — 0.005 (for the fermionic action 
with this small mass and corresponding to the 
same RGT as above), as a typical example. The 
solid line is the exact spectrum, the dashed is for 
the action truncated to the hypercube, the long- 
dashed is the truncated one with tiny corrections 
to satisfy the Symanzik condition, the dot-dashed 
is for the Wilson fermions. Fig. [g] shows the same 
quantities for to = 1.28. (Note that the Fermi- 
lab action [|l0| on the isotropic lattice gives a far 
worse dispersion relation: for m = it is complex 
and for ma > the kinetic mass is too small, see 
a figure in Ref. |l^. Remember, however, that it 
is assumed to be used on anisotropic lattices.) 

The spectrum of the full FP action is, of course, 
exact for any value of the fermion mass (or rather 
ma). The coefficients and A in eq. (|l|) depend 
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Figure 6. Fermion spectrum at ma ~ 1.28. 



on the mass, and they decay exponentiaUy in the 
distance r as oc exp(— 7r). A truncation at some 
small distance should necessarily induce 'ghosts' 
in the spectrum (higher branches, which are ar- 
tifacts) with an energy 7. A truncation is sen- 
sible therefore only if ma < 7. It is not yet clear 
what is the largest mass this condition allows (7 
will also depend on ma), but ma « 2 is still def- 
initely allowed. To study higher lattice masses, 
one might need to consider anisotropic lattices, 
anyhow. 

The mass dependence of the coefficients in 
eq. (|l^) follows a smooth curve which could 
be parametrized easily. In Fig. the ratios 



A(r, ma)/A(r, 0) and pi{r,ma)/ pi(r,0) are shown 
for different r values on the hypercube, only as 
an illustration to this statement. (Note that this 
dependence is obtained by an optimized 'running 
RGT'.) 
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Figure 7. Mass dependence of the ratios 
A(r, ma)/X{r, 0) and pi (r, ma)/pi{r, 0). 

As we saw, thermodynamic quantities are very 
sensitive to the cut-off effects. In Fig. p| for the 
case of massless fermions the ratio p/T is plot- 
ted vs. Nt, for the Wilson action, the full FP 
action and the FP action truncated to the unit 
hypercube (but still uncorrected for the tiny vio- 
lation of the Symanzik condition). Note that the 
full FP action has a correction vanishing expo- 
nentially with increasing Nt. 

9.3. Fermions and gauge fields 

The final goal of the improvement approach is 
to produce a good action of the full QCD. This 
goal is not achieved yet, but some important steps 
are already made. I shall review here the situa- 
tion within the RG approach. To derive the FP 
action for the case of fermions -I- gluons the fol- 
lowing steps have to be performed: 

1. Find the FP action for the pure gauge 
fields and free fermions, eqs. (|,|l|), and ob- 
tain the corresponding minimized fine fields 
from the coarse ones, U = U{V), ^ — V'(*), 
etc. 
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Figure 8. p/T'^ vs. Nt for massless fermions. 



The MIT group made a preliminary study |14[] 
in the quenched approximation — using a naive 
gauging of the free fermionic FP action, they mea- 
sured the pion dispersion relation with zero bare 
quark mass, at /? = 5. Instead of a light pion 
they obtained a heavy one {rriT^a « 3), but with 
a surprizingly good dispersion relation. 

Note that based on the 1-loop perfectness, for 
the full QCD FP action we do not expect a 
large perturbative renormalization. This should 
be, however, the consequence of terms non-linear 
in the link variables U in the interaction vertex 
RA{n,n',U) (cf. eq. (10)). Only after including 
these non-linear terms should one expect a small 
bare mass renormalization. 



2. Substitute this into the interaction part in 
the saddle point equation: 

5^(§,^',1^) = 5f(V5,^,C/)+ (20) 

rF(^,vi/,y,v;,^,c/)|^^^(^),^^^(^)... 

This equation is already satisfied for F = 1 
since it is the FP equation for the free fermionic 
action. For smooth fields V one can solve it by 
iteration, obtaining the linear (in Afj,{x)) approx- 
imation to the interaction vertex. This can be 
expressed in the gauge invariant form 

4>in)TARAin,n',U)^in') (21) 

where Ta is a Dirac matrix and RA{n,n' ,U) is 
a sum of weighted paths connecting the sites n 
and n' . For the MIT blocking this approximation 
has been done in Ref. Q. For the Bern/Boulder 
blocking(s) the details have not been published 
yet, but the main points are presented at this con- 
ference ||l2[. The interaction vertex i?^(n,n',[/) 
in the linear approximation looks roughly as fol- 
lows. For r = 1 or 7^ the coefficients are obtained 
by gauging the free fermion FP action: the near- 
est neighbours are connected by the direct link 
and the 6 staples with a weight of 50-50%, the 
others on the unit hypercube predominantly by 
the average of the corresponding shortest paths. 
The F = (7^1^ term (not present in the free fermion 
action) is rather small, the 75 and 7/^75 terms are 
even smaller. 



10. CONCLUSION 

The size of the lattice artifacts depends on the 
model investigated, on the particular form of the 
action used, on the quantity one measures, and, 
of course, on the actual lattice spacing a. The an- 
swer to the question that which improved action 
should be preferred depends also on many fac- 
tors — in particular the overhead in the MC sim- 
ulations, the expected improvement, etc. Since 
the computing costs grow like a (sometimes large) 
power of the inverse lattice spacing, even a large 
overhead can be tolerated if the corresponding ar- 
tifacts are sufficiently suppressed. 

At lattice spacings used in present numerical 
simulations the tree-level or one-loop Symanzik 
improved actions do not seem to improve suffi- 
ciently. One solution is to modify them in a non- 
perturbative way. Both the tadpole improvement 
and the method of non-perturbative determina- 
tion of the improvement coefficients go in this 
way. The actions used are relatively simple mod- 
ifications of the standard action, hence the over- 
head is relatively low. 

The more ambitious RG approach has a longer 
incubation period — it is not easy to derive and 
optimize the FP actions. Since their structure is 
more complicated the overhead is larger. The ex- 
perience collected so far indicates, however, that 
the resulting artifacts are very much suppressed, 
and it might compensate the extra complications 
of this approach. As mentioned, the trunca- 
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tion of the FP action can be made to preserve 
the Symanzik conditions, so it suggests a natural 
choice in the Symanzik improvement program. In 
measuring the topological susceptibility this ap- 
proach has the advantage to represent properly 
the classical instantons, avoiding the main prob- 
lems of the standard approach. One can specu- 
late that in cases where topological effects have 
strong influence on the behaviour of the system, 
the use of the FP actions is preferable in spite of 
the extra overhead. 

For spin systems, pure gauge theories and free 
fermions there exist good choices of improved ac- 
tions, which have passed various tests. For full 
QCD at present there are only promising candi- 
dates — in the coming year(s) the situation will 
be, hopefully, further improved. 

I would like to thank Anna Hasenfratz, Peter 
Hasenfratz and Peter Weisz for useful comments. 

This work was supported by the Swiss National 
Science Foundation. 
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